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ABSTRACT 

We revise the physical conditions in the O"'""'" zone of the planetary nebula 
NGC 6543, obtaining two different estimates of the electron temperature (Tg) 
and one estimate of the electron density (iVg). The electron temperature is com- 
puted by means of two independent methods, the nebular-to-auroral ratio [O ill] 
A5007/A4363, and the diagnostic diagram that combines A 5007 to the [O ill] 
infrared lines 52 /z and 88 /x. The optical and infrared fluxes have been obtained 
from archive HST/WFPC2 images and ISO LWS spectra respectively, and the 
continuum intensity in the optical has been measured from narrow-slit spectra 
obtained with the Isaac Newton Telescope at La Palma. The measured con- 
tinuum intensity is higher than predicted by recombination theory under the 
hypothesis that all the Lya photons either escape or are destroyed. This fact 
can be explained in terms of an enhancement of the 2-photon continuum due to 
Ly-a conversion, a process that depends strongly on the local structure of the 
nebula. Alternative possibilities, outside the framework of recombination theory, 
have also been considered: e.g., the optical tail in the X-ray emission of a very 
hot plasma, and dust scattering of stellar radiation, but these hypotheses are not 
supported by quantitative estimates. 

While the electron temperature and density derived from the diagnostic di- 
agram agree with the most recent determination, the temperature derived from 
A 5007/ A 4363, Te(0++)opt, is somewhat smaller than previously published values. 
We discuss several technical issues that contribute to the overall uncertainty in 
our results, focusing on the instrumental effects that might bias the [O ill] A 4363 
intensity. We also discuss the effects of the coUisional de-excitation of the 0+"*" 
ground terms on the relation between Te{0^^)opt and Te{0^^)iR. 
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1. Introduction 

Knowledge of physical conditions {Tg and Ng) is an essential ingredient in the under- 
standing of planetary nebulae. Information on the structure, chemistry and kinematics of 

these objects comes from individual diagnostics and from the cross-comparison of different 
diagnostics. In this work, two different diagnostics will be used to investigate the physi- 
cal conditions in the planetary nebula NGC 6543: the temperature-sensitive [O III] ratio 
J(A5007)//(A4363), and the (T^, N^) diagnostic diagram based on the [O III] optical and 
infrared lines, which was first proposed by Dincrstcin, Lester, & Werner (1985) (henceforth 
DLW; see also Dinerstein 1983). 

The first method, the one based on the /(A 5007)//(A4363) ratio, is one of the most 
popular ones. Advantages of this method are that the lines involved are bright, and close 
enough to be easily observed with the same instrumental setup; a caution to be taken is that 
this ratio selectively samples the hottest zones, so that in presence of strong temperature 
fluctuations the derived temperature, Te(0''"+)opt, is biased towards the highest values (see, 
e.g., Peimbert 1995). Furthermore, this ratio measures the temperature of the 0+''" region, 
which is not necessarily representative of the whole observed volume. The second method, the 
diagnostic diagram by DLW, allows the simultaneous determination of the electron density 
and temperature by combining the optical and the FIR fine-structure [O III] lines. The 
two infrared lines [O III] 52fi and 88fi arise from low-energy levels of the ground term, so 
their ratio 52/i,/88/i is almost insensitive to the temperature, while the critical densities of 
their upper levels are different, so the ratio is strongly sensitive to density. On the other 
hand, a ratio like A 5007/52// is strongly dependent on temperature, due to the very different 
excitation energies of the upper levels of the two lines, and it has a weak dependence on 
density too, due to the difference in the critical densities. By combining the two line ratios, 
one can obtain a simultaneous measurement of the electron temperature, Te{0~^~^)jji, and 



^Based on observations with ISO, an ESA project with instrument funded by ESA Member States and 
with the participation of ISAS and NASA. 
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density, Ne{0'^'^)iR (we will use throughout the text the subindex IR when referring to these 
Tg and determinations, though an optical line is involved also; this choice responds to 
simplicity criteria, and highlights the fact that the density is, to a good approximation, a 
function of the ratio of the two infrared lines alone). One of the main advantages of this 
method is its internal consistency, since all the lines involved belong to the same ion, so 
that the same volume is mapped to derive both temperature and density (although the 
optical and infrared emissivities might follow different distributions within the O"^"*" zone if 
temperature or density fluctuations are present). Two major disadvantages are consequences 
of the lines involved lying far apart in the spectrum: first, the lines must be observed with 
different instruments, a procedure implying the risk of systematic offsets in the calibration; 
and second, any uncertainty in the reddening correction is amplified by the large wavelength 
baseline. 

By means of these methods, we want to reassess the physical conditions in the bright core 
of NGC 6543. A similar analysis was performed on a small sample of PNe, including NGC 
6543 itself, by DLW. These authors used infrared data obtained with the Kuiper Airborne 
Observatory (KAO), and complemented them with optical spectra from the Lick Observatory 
Crossley reflector; the analysis was later repeated by Dinerstein, Haas, Erickson, & Werner 
(1995) (henceforth DHEW), with quite different results. In the work we present here, the 
intensity of the infrared hnes is obtained from ISO archive spectra, and the intensity of 
[O III] A 5007 and [O lll] A 4363 from HST archive images. Additionally, Ha and H/3 HST 
images were used to determine the reddening correction to be applied to the other images, 
and narrow band HST images around A 6584 were used to correct the Ha images for the 
[N ll] contribution. Finally, long-slit, ground-based data obtained with the Isaac Newton 
Telescope at the Obscrvatorio del Roque de los Muchachos were used to correct the HST 
images for the contribution of other lines and the continuum. These data will be described 
in Sections 2, 3 and 4. 



1.1. NGC 6543 

NGC 6543 (the 'Cat's eye') is a very well known, low-excitation planetary nebula. Many 
studies have been devoted to understand its complex structure; an incomplete list of the most 
recent ones includes Reed et al. (1999); Balick, Wilson, & Hajian (2001); Hyung et al. (2001). 
The bright core of the nebula can be roughly described as the superposition of two ellipses, 
equal in shape but perpendicular to each other (Miranda & Solf 1992). This region is about 
25" in diameter, and can be seen in direct images in the Balmer and [O ill] lines. Other 
structures form part of the nebula: surrounding the bright core, a series of dim concentric 
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rings is visible in [O ill], extending at least out to 70" (Balick, Wilson, & Hajian 2001); in 
turn, the rings are surrounded by a giant filamentary halo of about 300" in diameter (Balick, 
Wilson, & Hajian 2001; Middlemass, Clegg, & Walsh 1989); X-ray emission coming from 
the bubble surrounding the central star has been recently reported by Chu et al. (2001); the 
variability of the central star has been investigated by Balick, Wilson, & Hajian (2001). 

2. The ISO data 

NGC 6543 was routinely observed by ISO for calibration purposes, thereby a large set 
of spectra is available. We selected 20 out of the 92 Long- Wavelength Spectrometer (LWS) 
spectra of NGC 6543 available in the ISO data archive, all of them in the AOT LOl observing 
mode; the selected spectra were taken during two separate periods, close to the beginning 
and the end of ISO operational life respectively (see Table 1). A detailed explanation of the 
LWS structure and data sets is provided in The ISO Handbook, vol. IV, LWS - The Long 
Wavelength Spectrometer (Gry et al. 2001). 

Each spectrum is comprised of 10 independent wavelength ranges, which correspond 
to the 10 detectors that form the LWS. These wavelength ranges overlap partially. Both 
emission lines considered by us fall in the overlap region of two consecutive ranges: i.e., 
[O III] 52/x falls on both the SWl and SW2 detectors, while [O ill] 88/x falls on both the 
SW5 and LWl detectors. Additionally, the spectrum of each of these 10 wavelength ranges 
is composed of 6 sub-exposures, that correspond to independent scans of the detector. Thus, 
we can in principle perform, in our sample, (20 spectra) x (2 detectors) x (6 scans) = 240 
independent measures for any of the two emission lines of interest, or 120 with each detector. 
In practice, the number of independent measurements of a given line with a given detector 
is smaller than 120, due to the following reasons: in a few cases, one or two out of the six 
scans of a given exposure do not span the whole detector range, failing to include the line; in 
other cases the quality of the spectrum is very poor, so that the line-intensity measurement 
is unfeasible. However, the measurements are still numerous enough (77 on average per line 
and detector) to adequately sample the data distribution. 

The data were reduced using the latest version of the ISO automatic analysis pipehne 
available in November 2001. The flux measurements were performed by model gaussian 
fitting. Figure 1 represents the continuum flux versus emission line intensity for all the 
measurements of the two features, 52yU and 88/i. In the upper plot, the filled dots and the 
open circles represent the data taken with the SWl and the SW2 detectors respectively; in 
the lower plot, the dots and the open circles represent the data taken with the SW5 and 
the LWl detectors respectively. The ellipses delimit the region within 1-a from the median 
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values. 

Table 2 summarizes the median and dispersion values of the continuum and line mea- 
surements for the four detectors. The data of each epoch are shown both as individual data 
sets, and combined together. The following conclusions can be drawn from the analysis of 
the data sample: 

a) the SWl and SW5 continuum data are much more spread out than the corresponding 
SW2 and LWl data; 

b) the median intensity values of the 88// line obtained with SW5 and LWl are in excellent 
agreement, whereas the median intensity values of the 52// line obtained with SWl and 
SW2 are more discrepant, although they still agree at the 1-a level; 

c) the median intensity values of the 52/i line obtained with SWl suffered a drift of more 
than la from Epoch I (June 1996 + July 1996) to Epoch II (December 1997 + January 
1998). This long term drift is known and documented in the ISO LWS Handbook. 

Based on these results, we decided to perform our analysis on the hue intensities mea- 
sured with the SW2 and LWl detectors only, neglecting the data obtained with the other 
two detectors. The exclusion of SWl is primarily motivated by the drift in the hue flux data, 
which make them unreliable. On the other hand, the exclusion of SW5 is only motivated by 
the larger spread in the continuum data with respect to those of the LWl detector, a feature 
that does not necessarily imply a lower quality of the line fluxes obtained with this detector. 
In fact, since SW5 has a spectral resolution twice that of LWl, the recorded dust continuum 
level per resolution clement is lower, making this detector more sensitive to the dark current 
subtraction: as a result, a fraction of the larger spread in the SW5 continuum data may be 
caused by the uncertainties in the subtraction of the dark current, an uncertainty not affect- 
ing the line fluxes. Nonetheless, we preferred to adopt a conservative approach and exclude 
both SWl and SW5 from our analysis. With such a choice, the adopted line intensities are: 

/(52//) = (5.10 ± 0.61) X 10"^° ergsec"^cm-2 (1) 

and 

/(88/f) = (1.39 ± 0.17) X 10-^° ergsec-^cm-2, (2) 

where the quoted value is the median of the sample, and the uncertainty is the dispersion 
in the data, which does not include possible systematic error sources, such as, e.g., the flux 
calibration. These values are in excellent agreement with those reported by Liu et al. (2001), 
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/(52/i) = (5.11 ± 0.13) X 10-^° erg sec"^ cm-^ and 7(88//) = (1.49 ± 0.04) x 10"^° erg sec"^ 
cm~^, although Liu et al. (2001) quote rather small uncertainties. 

A flux correction factor should be applied to take into account the extension of the 
source, since the ISO LWS analysis pipehne is calibrated on Uranus and it assumes that 
the target is a point source. This flux correction factor accounts for the relation between 
the ISO telescope point spread function and the LWS instrumental profile. DLW applied 
the corresponding flux correction to their KAO data, a correction that was quite small, 
amounting to only 4% for the 52 /i line. In our case, the fraction of flux diffracted out of the 
aperture is even smaller, and it is therefore negligible with respect to the other uncertainties 
involved in the process, such as, e.g., the flux calibration itself. Consequently, we did not 
apply any correction for this effect. 



3. The Isaac Newton Telescope long-sIit data 

NGC 6543 was observed spectroscopically as part of a wider study of PNe on 1995 July 8 
and 9. We used the Intermediate Dispersion Spectrograph attached to the 2.5m Isaac Newton 
Telescope (INT) at the Observatorio del Roque de los Muchachos, on La Palma. The 4' x 1" 
slit was placed at position angles 5° and 20°, and included the central star. The R1200B 
grating and the TEK3 CCD in the 235 mm camera provide a dispersion of 0.85 A/pixel 
and a spatial sampling of 0.71"/pixel. The seeing was stable at 1.1", and the integration 
time was 1000 s. The spectra were calibrated using the standard steps of ffat fielding, bias 
subtraction, wavelength calibration, and flux calibration using flux spectroscopic standard 
stars observed during the same night and with the same setup, but through a 10 arcsec 
wide slit; the spectra were also corrected for light dispersion inside the spectrograph. The 
inset in the top panel of Figure 2 illustrates the very high quality of these data, showing a 
logarithmic plot of a sample nebular spectrum in the wavelength range 3400-4500 A. These 
data form part of a much larger data set that will be presented elsewhere. 



4. The HST images 

The HST data archive at ST-ECF was queried for WFPC2 images of NGC 6543. We 
retrieved images through the narrow band fllters F437N, F487N, F502N, F656N, and F658N, 
centered on the emission lines [O ill] A 4363, H/5, [O ill] A 5007, Ha, and [N ll] A 6584 re- 
spectively. Table 3 summarizes these data sets. For each filter, the different exposures 
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were combined with the IRAF^/STSDAS command crrej, to obtain a final image at each 
waveband free of cosmic rays. In each of these images, the fluxes were measured within the 
PC frame of the WFPC2. Similar analyses have been performed by Lame, Harrington, & 
Borkowski (1997) and Hyung et al. (2001) for NGC 6543, and Rubin et al. (2002) for NGC 
7009. Since our results differ from those of both Lame, Harrington, & Borkowski (1997) 
and Hyung et al. (2001), we will describe our procedure in full detail, so as to allow fu- 
ture identification of the causes behind the discrepancies. The calibration is conceptually 
straightforward, but the procedure contains several potential pitfalls, such as, e.g., the ze- 
ropoint calibration, the adoption of the proper width for the filter, or the estimation of the 
contribution to the fiux of the continuum and of neighbouring lines. 



4.1. Photometric zeropoint calibration 

The zeropoint of an instrument is defined as the intrinsic fiux of an object producing 
one count per second. The raw data number (DN: the ratio between the electron count in 
each pixel, and the gain factor, which is 7 in our case) obtained from the HST images must 
be multiplied by the appropriate zeropoint in order to obtain the flux in energy units. The 
zeropoint for the WFPC2 is codified in the header word PHOTFLAM. As many other cali- 
bration parameters, PHOTFLAM has been periodically modified, as the calibration routine 
improved. The PHOTFLAM values can be found either in the header of the images, or 
in the HST-WFPC2 Data Handbook (Baggett et al. 2002). Unfortunately, the two sources 
give different values; the difference is quite sizeable (almost 4%) in the case of the Ha filter 
F656N, and smaller than 1% for the other filters. As the HST helpdesk advised us to do, 
we sticked to the PHOTFLAM values of the WFPC2 Data Handbook, which are fisted in 
Table 3. 



4.2. Width of the filter 

The narrow WFPC2 filters are characterized by several parameters. One of them is the 
filter's rectangular width RECTW, defined as the width in A of a rectangle with the same 
total area as the total transmission curve, which describes the transmittance of the system, 
and height given by the peak in the curve. Since the DN is given per unit wavelength, the 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 



-8- 



data measured in the WFPC2 images must be multiplied by the RECTW value of the cor- 
responding filter. As already noted by Rubin et al. (2002), RECTW should not be confused 
with the SX values given in the WFPC2 Instrument Handbook (Biretta et al. 2001, Appendix 
1). RECTW can be computed for each filter with the IRAF task synphot/bandpas. The 
values for the filters considered in this work are given in Table 3. 

4.3. Field of view 

The PC frame is a square of 800 x 800 pixels, corresponding to about 37" x 37". 
The central 40 x 38 pixels were excluded, to avoid contamination by the central star. The 
outermost pixel rows and columns arc affected by optical problems related to the spherical 
aberration of the telescope. Due to the low photometric quality of these pixel rows they 
were excluded from our analysis, and the fiux was measured in a square of 690 x 670 pixels, 
corresponding to about 32" x 30". This square is smaller than the field of view of the ISO 
LWS, which has a diameter of about 80" at the wavelengths of interest. However, the flux 
emitted by the nebula is completely dominated by the innermost 30", which accounts for 
more than 97 percent of the total (Figure 3), and the two instruments effectively sample the 
same area (this is also shown, e.g., by the brightness distribution of the [Ne III] line at 15 ^: 
see Pcrsi ct al. 1999). This issue is also addressed by DLW, who match the KAO 40 — 50" 
aperture by taking the ground based observations of [O IIl]A5007 and B.(3 through 45" and 
70" diameter apertures, and using the [O ill] A 4363 fluxes by Bohuski, Dufour, & Osterbrock 
(1974) through a similarly wide aperture. 

4.4. Flux contribution of neighbouring lines 

A fraction of the flux measured in each fllter is contributed by lines other than the 
principal line. In the case of the four fllters we considered, F487N and F502N contain only 
R/3 and A 5007 respectively whereas F437N includes H7 A 4340 and He I A 4388, and F656N 
includes [N ll] AA 6548, 6584. In this section, we describe the procedure followed to subtract 
these contributions from the total measured fluxes. 

The total flux of [N II] A 6584 was computed from an image taken with the F658N fllter, 
and then subtracted from the F656N image taking into account the transmission efficiency of 
the filter at the position of A 6584. The [N ll] A 6548 contribution was computed theoretically 
from the fiux of [N ll] A 6584. Since the F658N filter is, in turn, contaminated by Ha, the 
decontamination procedure was iterated to give decontaminated Ha and [N ll] A 6584 fiux 
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values. After this procedure was applied, the total Ha flux lowered by about 3%. 

The total flux in the F437N filter includes contributions of H7 A 4340 and He I A 4388. 
Both contributions are difficult to estimate, since they depend on the exact position of the line 
in the filter. In the case of He I A 4388, the intensity of the fine itself is also unknown, whereas 
the total intensity of H7 can be computed from I (Ha). We estimated these contributions 
multiplying the nebular spectrum of the long-slit spectra (with the central star excluded) 
by the F437N transmission curve, and then measuring separately each component on the 
resulting pattern. The top panel in Figure 2 shows the nebular spectrum of NGC 6543 along 
position angle 5° in the wavelength region around H7 4340, [O III] 4363 and He I 4388. The 
dashed line is the actual spectrum, while the full line is this spectrum convolved with the 
HST F437N filter response (the total system peak response of 0.03 (Biretta et al. 2001) has 
been divided out so as to make the two spectra comparable). The inset shows a sample of 
the spectrum along the sht in a logarithmic scale, to provide an indication of the high S/N 
of this data set. We found that H7 and He I A 4388 contribute on average 12% and 4% 
respectively to the total flux in F437N, while [O ill] A 4363 provides on average 26% of the 
flux, with a typical deviation of 3%. This is plotted in the bottom panel of Figure 2, which 
shows the contribution of A 4363 to the total intensity in F437N, normalized to 100 (dotted 
circles, left handside scale; the central ±3" are excluded from the analysis); the rest of the 
fiux detected through the filter is contributed by H7 and He I 4388, other very faint emission 
lines, and by the underlying nebular continuum (see top panel). Also plotted for reference is 
the variation of the [O ill] 4363 line intensity along the slit (solid fine, right hand side scale) . 



4.5. Continuum subtraction 

The underlying continuum contributing within each of the HST images can be ascer- 
tained either theoretically or empirically. The first method requires a detailed computation 
of the efficiency of all the processes involved, while the second consists in measuring on a 
spectrum the relative importance of each line with respect to the underlying continuum, 
with the filter's transmittance curve taken into account. We chose to do it both ways, since 
the comparison between the results illustrates quantitatively the difficulties of an accurate 
theoretical computation of the nebular processes. 
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4-5.1. Theoretical continuum contribution 

The processes contributing to the nebular continuum flux are H I, He I, and He II recom- 
bination, bremsstrahlung, and 2-photon decay. The relative importance of these processes 
depends on both the wavelength and the nebular electron temperature, and on the He I 
and He ll abundances. Table 4 shows the frequency dependence of the continuum-emission 
coefficient, 7,^, for each process at the three wavelengths of interest, in units lO^^^ erg cm^ 
sec~^ Hz~^; the data have been obtained interpolating in the {u, log j,^) plane between the 
values hsted by Brown & Mathews (1970) at Tg = 8000 K. 

The 2-photon continuum intensity is quite difficult to compute accurately, as it depends 
not only on the local physical conditions Ne and Te, but also on the fate of the Lya photons 
produced by recombination. In optically thin or dusty nebulae, all the electrons on the 2^P 
level produce a Lya photon that eventually escapes from the nebula or is destroyed by dust, 
whereas those on the 2^S level decay to the ground state producing 2-photon continuum, 
the population on 2^S and 2^P levels being determined by the recombination cascade and 
the mutual collisional transitions. In practice, Lya photons can be scattered repeatedly 
before they escape from the nebula, and at each scattering the absorbing electron has a 
finite probability of suffering a 2^P — > 2^S transition and decaying via 2-photon emission. 
Thus, with given local physical conditions {Ng, Tg), the actual 2-photon emission can vary 
from a minimum value (corresponding to the case in which all the Lyct photons are either 
destroyed by dust or escape from the nebula) to a maximum value (corresponding to the 
case in which all the Lya photons are eventually converted to 2-photon continuum). In 
Table 4 and Figure 4 we represent, as a function of the wavelength, the two limiting values 
for the 2-photon continuum, and the corresponding values for the total continuum. The 
impact of the occurrence of Lya conversion on the total continuum level depends on the 
wavelength considered; in the case of the wavelength range around A 4363, where 2-photon 
emission is a dominant process, Lya conversion can enhance the total continuum by as 
much as a factor 2.5. This would be the asymptotic case of a very thick, dust-free nebula, 
and represents just an upper limit to be considered in the following. We assumed in the 
calculations Tg = 8000 K and Ng = 1650 cm~^; the density corresponds to the observed 
52yu/88/i ratio: see Section 5. The collisional rates listed by Osterbrock (1989) for the 2^S 
— 2^P transition have been linearly extrapolated in temperature down to Tg = 8000 K; this 
might seem a rough approximation, but our results are fairly insensitive to it, both because 
the collisional data at = 10000 and 20000 K are quite similar, and because collisions are 
outnumbered by spontaneous decays at these nebular conditions anyway. 

The He II recombination continuum is not listed because it is not relevant for NGC 
6543, due to the low He++ abundance (A^(He++)/A^(H+) ~ 10-^ Hyung et al. 2000). On 
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the other hand, A'"(He"'")/A^(H+) ~ 0.11, so that the He I recombination continuum cannot 
be neglected. 

The continuum-emission coefficient j^,, measured in erg sec~^ cm~^ Hz~^ ster~^, is re- 
lated to through ji, — A^+A^e7i//47r. Figure 4 shows ji, normalized to A^(H+)A^e/47r, i.e., the 
He I recombination contribution has been multiplied by the factor A'"(H+)/A^(He"'"), so that 
the contributions of the three processes to the continuum can be directly compared. The 
specific intensity emitted through each process depends on the width of the filter's passband 
RECTW: 4 = /rectw>^^ = ^+^e/47r /j^g^^^ 7^(ii/. 

The continuum contribution is very important in filters centered on weak lines, such 
as [O III] A 4363, but it barely affects the intensity in filters centered on strong lines, such 
as Ha or [O ill] A 5007. Our computations show that, at Tg = 8000 K and neglecting 
Lya conversion, the continuum yields 40% of the total intensity in the F437N filter around 
[O III] A 4363, 0.3% of the total intensity in the F502N filter around [O ill] A 5007, and 
0.5% of the total intensity in the F656N filter centered on Ha. As anticipated before, 
we also measured the continuum contribution on long-slit spectra, and compared it to our 
theoretical computations. Our findings are discussed in next section. 

4-5.2. Continuum contribution from direct measurement 

The expected continuum intensity can be compared to direct measurements of the con- 
tinuum level in an optical spectrum. We followed this procedure measuring (through the 
corresponding filter responses, as described in section 4.4) the integrated continuum flux and 
relevant fine intensities in the two long-sht spectra of NGC 6543 described in Section 3. By 
means of this procedure, we find that the continuum contributes 58% of the total intensity 
in the F437N filter around [O ill] A 4363 (cf. Figure 2), significantly more than expected on 
the basis of our previous calculations. 

Three likely candidate mechanisms that would enhance the continuum level are X-ray 
emission, dust scattering, and 2-photon emission. X-ray diffuse emission has been recently 
reported from high resolution spectral images by Chandra (Chu et al. 2001); this is produced 
in the mass-loading shocked region where the fast stellar wind interacts with the inner 
gaseous shell. The optical tail of such emission could be the origin of the observed extra 
continuum. To test this point, we modeled the central cavity of NGC 6543 by means of 
simple, spherical hot-plasma models at different temperatures, normalized to the soft X-ray 
flux observed by Chandra, and computed for each model the contribution to the F437N 
band. Large uncertainties in this calculation are introduced by the values assumed for the 
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plasma temperature, the exact lower energy limit seen by Chandra (which is around 0.1 - 
0.2 keV), and the intervening absorption by neutral gas (which we did not take into account 
in our estimates). For example, a T = 1.7 x 10^ K plasma, normalized to the emission 
in the (0.1 - 10 keV) interval, produces only 0.001% of the observed extra-continuum in 
F437N, while a T = 5 x 10^ K plasma, normalized to the emission in the (0.4 - 10 kcV) 
interval, produces 10% of the observed extra-continuum in F437N. Although these are only 
back-of-the-envelope computations, altogether they suggest that the extra continuum cannot 
be explained in terms of hot plasma emission. 

A different mechanism that could explain the extra continuum is dust scattering of the 
central star's radiation. Even if the external screen extinction is rather small, C(H/3) = 0.07 
(see below) , a cursory inspection of SWS-I-LWS archival spectra of NGC 6543 shows a thermal 
bump due to hot dust peaking at ~ 30//. Although this mechanism is a plausible candidate 
for the origin of the extra continuum, the detailed modeling of the internal dust distribution 
in NGC 6543 is beyond the scope of this paper. 

A third possibility is 2-photon emission with enhanced efficiency. As discussed in Sec- 
tion 4.5.1, in our estimate of the 2-photon emissivity we made use of a simplified treatment 
(Osterbrock 1989; Brown & Mathews 1970), which neglects Lya conversion. The efficiency 
of Lya conversion to produce 2-photon emission varies between in optically thin or dusty 
nebulae (all the Lya photons either escape or are destroyed by dust) and a maximum (all 
the Lya photons are eventually converted to 2-photon emission). The detailed calculation of 
the actual efficiency of this process is virtually impossible, since it would require knowledge 
of the local escape probabilities and the dust structure of the nebula, and the implementa- 
tion of this information in a 3-D computation. The only feasible alternative is to bracket 
this process by computing the 2-photon emission in the extreme case of maximum efficiency. 
Following the formalism of Brown & Mathews (1970), it is found that the 2-photon emission 
at A 4363 can be enhanced by Lya destruction by up to a factor of 3.5, rising the total con- 
tinuum by up to a factor of 2.5. Figure 5 compares the observational data to the continuum 
range predicted by recombination theory, as a function of the H7 intensity. The data are 
compatible with a Lya-conversion enhancement of the 2-photon continuum, and suggests 
that it is not necessary to invoke other processes to explain the measured continuum level 
(the only discrepant points are some of those with the lowest intensity values, i.e. with the 
largest errors). The figure also shows that the relation between H7 and the measured con- 
tinuum is not linear, and that there are large fluctuations around the best-flt straight line, 
as expected if the continuum in excess of the minimum level arises from Lya conversion: 
stated otherwise, the Lya conversion process acts with different efficiency across the nebula. 
This figure demonstrates that it is not possible to rely on theoretical calculations alone to 
correct the total flux in each image for the continuum contribution: flrst because we would 
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ignore, a priori, the average efficiency of tfie Lya conversion process (tlie slope on tfie H7- 
continuum relation, which is provided by the slit spectrum); second because, even with this 
piece of information available, wc would still ignore the large fluctuations of the continuum 
around such average level. These fiuctations naturally arise as a consequence of the local 
character of the 2-photon enhancement process, and, since theory alone provides no way 
of predicting where and how large they are across the image, they would translate directly 
into an uncertainty in the measured A 4363 flux. Furthermore, this uncertainty would be 
relatively much larger for A 4363 than it is for the continuum, since the relative contribution 
of the hne in the F437N filter is less than one half that of the continuum, so that a typical 
uncertainty of, e.g., 20% in the continuum would translate into an uncertainty of 45% in the 
line. This would be much larger than the typical uncertainty that can be read directly from 
Figure 2, and it would not avoid the need for a spectroscopic measurement; so this strategy 
to determine /(A 4363) across the image is definitely not viable. 

Note also that this conclusion holds even under the hypothesis in which the 2-photon 
continuum is at its minimum level (the only case that would allow to predict the contin- 
uum produced by recombination by means of a recombination line) and that the additional 
continuum flux comes instead from a different source: as a matter of fact, the continuum 
is not proportional to H7, hence it cannot be predicted theoretically, no matter which pro- 
cess is responsible for the departure from direct proportionality. An important consequence 
is that no temperature maps can be obtained using this kind of data, unless very specific 
circumstances occur: one would be that A 4363 dominate the light in the filter, so that only 
a rough correction for the continuum would suffice; a second one would be the availability 
of tunable filters to measure directly the continuum level near A 4363, bypassing the need 
for theoretical predictions. The lesson to be learned from these computations is that pre- 
cision photometry of weak lines cannot be done with the HST/WFPC2 in the absence of 
data specifically designed to measure the continuum. As a matter of fact, an additional 
uncertainty is introduced even by the estimation of the contribution of the various lines 
inside each filter, unless a reliable calibration is available. Both these problems and possi- 
ble solutions to them - continuum and line subtraction - have been addressed by O'Dell & 
Doi (1999) and O'Dell, Peimbert, and Peimbert (2003); their arguments will be described 
in Section 6.1. To proceed with the discussion, we will therefore assume that the A 4363 
percentage within the filter F437N is the one provided by the bottom panel of Figure 2, i.e. 
/(A 4363)//tot(F437N) = 0.26 ± 0.03. It must be emphasized, however, that in the absence 
of slit spectra or sophisticated theoretical tools to estimate a priori the continuum intensity, 
such uncertainty would be much larger, amounting to the difference between the minimum 
and maximum continuum value admitted by recombination theory. 

In their paper on NGC 6543, Hyung et al. (2001) pointed out inconsistencies between 
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the theoretical estimation and the continuum level derived from observations (in their case, 
ground-based cchelle spectra). It is not clear, however, whether such inconsistencies can 
be tracked down to what we point out in this work, as their procedure differs from ours 
in at least two relevant aspects: first, the theoretical continuum flux computed by Hyung 
et al. (2001) includes only the contributions by H I recombination and bremsstrahlung. 
However, Fig. 1 in Brown & Mathews (1970) shows that, at Tg = 10000 K, the continuum is 
dominated by 2-photon continuum at A 4363, and this process is only slightly less important 
than H I recombination at A 5007. At — 8000 K the 2-photon contribution is even more 
important. The 2-photon continuum must therefore be included even in a rough estimation 
of the theoretical continuum flux. Second, Hyung et al. (2001) adopt smaller values than 
we did for the width of the filter: 25.2 A, 26.9 A, and 21.9 A for F437N, F502N, and 
F656N respectively, instead of 31.829 A, 35.781 A, and 28.336 A. Apparently, these values 
correspond to the effective width of the bandpass 5\ defined by Biretta et al. (2001) (pag. 
149 and Appendix 1; see also Section 4.2 above); however, we might be misunderstanding 
their procedure, since this point is not explicitly stated in their paper, and their Ha width 
slightly differs from the one given by Biretta et al. (2001). But if we are understanding 
their procedure, then the difference in the adopted filter widths propagates directly into 
differences in the estimated continuum flux. This introduces a large bias in all the derived 
fluxes. Furthermore, since the ratios of different RECTW values are not equal to the ratios 
of different 5As, the intensity ratios of different images are correspondingly affected (the 
way it affects the derived A 4363 is more subtle, because the continuum contribution to be 
subtracted from the F437N flux is directly proportional to the adopted bandwidth). 

It seems probable that Lame, Harrington, & Borkowski (1997) also relied on a theo- 
retical computation to estimate the theoretical continuum in the F437N image, since the 
number they quote, Fc(F437N) = 0.015 x F(H/3), is comparable (actually, smaller) than our 
theoretical computation. However, their work is just a short conference contribution, and 
they do not specify how such estimate is obtained. 

A similar procedure was also followed by Rubin et al. (2002) in their analysis of the plan- 
etary nebula NGC 7009. These authors obtain a temperature map of NGC 7009 by means of 
WFPC2 images taken in the same filters considered by us. To estimate the line contribution 
inside the WFPC2 filters, they make use of HST/STIS data, modified to take into account 
the transmittance of the corresponding WFPC2 filters. Following this procedure, they find 
that the continua measured from STIS are higher than expected at the temperature of the 
nebula. They interpret such discrepancies as an artefact due to insufficient signal-to-noise in 
the STIS data, and therefore adopt in their analysis a continuum flux from theoretical com- 
putations; it is possible, however, that the problem they found with the continuum in NGC 
7009 is analogous to what is found in NGC 6543, namely 2-photon continuum enhancement. 
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In any case, the uncertainty on the continuum is less crucial in their work than it is in the 
case of NGC 6543, because A 4363 is relatively more intense in NGC 7009: they measure 
in F437N a A 4363/continuum ratio equal to 3.7, i.e. almost one order of magnitude higher 
than the one observed in NGC 6543. 

Summarizing this section on continuum, it is important to stress that there is no safe 
way to predict accurately the A 4363 contribution to the total flux in the filter F437N, unless 
this line is very bright. Because of this difficulty, special care should be given to this point, 
and the relevant assumptions should always be made explicit. 

4.6. Reddening correction 

The mean value of the Hq;/H/3 flux ratio across the HST image of NGC 6543, with 

the [N II] and the continuum contributions subtracted, is Ha/H/?=3.06±0.12 (cf. Figure 6). 
Adopting the extinction law by Cardelli, Clayton, & Mathis (1989) with /?v/=3.0, and fitting 
to the theoretical Balmer decrement at Tg = 8000 K (Storey & Hummer 1995), a reddening 
coefficient C(H/3) = 0.07 ± 0.06 is obtained. 

The three C(H/3) values compiled by Hyung et al. (2001), C(H/5)=0.18, 0.20 and 0.30, 
are significantly higher than our result. We ignore the origin of this discrepancy, since none 
of the papers referenced by Hyung et al. (2001) lists the observed Ha fiuxes, nor do they 
provide details on the assumptions underlying the computation of C(H/3). On the other 
hand, our adopted C(Rj3) agrees with the value determined by other authors. For example, 
Kaler (1983) found C(H/3)=0.10, and the average value from the literature compilation by 
Kaler, Shaw, & Browning (1997) is (C(H/3)) = 0.09. 

5. Density and temperature in the O"'""'" zone 

The line intensities obtained by the measured fiuxes after dereddening and empirically 
correcting for continuum contribution, are: 

/(A 4363) = (3.95 ± 0.84) x 10"^^ ergsec"^ cm'^ (3) 

and 



/(A 5007) = (1.93 ±0.29) x lO"*^ ergsec"^ cm'^ 



(4) 
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The quoted uncertainties have been estimated as follows: the uncertainty on C(H/3), 
0.06 dex, translates into uncertainties of 0.07 dex and 0.06 dex for A 4363 and A 5007 re- 
spectively; the continuum subtraction adds another 0.05 dex to the uncertainty on /(A 4363) 
(Section 4.5.2); finally, further uncertainties should be added to account for the shift in 
wavelength of the bandpass with temperature, and the photometric errors (cf. Biretta et al. 
2001, Table 8.11); however, this is a very complex issue, as it will be shown in Section 6.1.1, 
and this source of uncertainty will not be considered for the moment. 

These values yield A5007/A4363 = 490±130, and Te(0++)opt = 7500 ± 450 K, a tem- 
perature lower than previously published values. For example, Hyung et al. (2001) find 
Te = 8000 ~ 8300 K, Kingsburgh, Lopez, & Peimbert (1996) find = 7950 ± 100 K, Peim- 
bert, Torres-Peimbert, & Luridiana (1995) find = 8334 K, and DLW adopt A 5007/ A 4363 
= 363 from Bohuski, Dufour, & Osterbrock (1974), which yields, using the same atomic data 
as us, Te(O++)opt=7900 K. 

Combining the infrared intensities measured on ISO spectra with the [O ill] A 5007 
intensity derived from the HST data, it is possible to simultaneously derive the temperature 
and density in the 0++ zone of NGC 6543 by means of the infrared-line diagnostic described 
by DLW. The diagnostic diagram is shown in Figure 7, where the point representing the 
data has been plotted. The data points by DLW and DHEW are also plotted. The values 
obtained are = 1650^4qq cm~^, Tg = 8600 ±500 K, while the corresponding values quoted 
by DLW and DHEW arc N^= lOOOOt^oo cm-^, = 5800 ±300 K and iVe= 2000t|[j|j cm'^, 
Te = 8500 ± 500 K respectively 

Figure 7 and Table 5 show that the point by DLW and the one we obtain here disagree 
by much more than la, both in Tg and Ng. The three lines are discrepant, particularly 52 fi 
and A 5007 (however, the discrepancy in A 5007 depends solely on the adopted C(H/9), since 
the difference between the measured fiuxes is less than 1%). These differences go in the 
same sense, so that the A5007/52/i ratio is not very different in the two cases: however, the 
upward bending of the isotemperature contours in the high-density region of the diagram 
enhances the difference between the corresponding temperatures. On the other hand, our 
results are in excellent agreement with those by DHEW, though unfortunately the fluxes of 
individual lines cannot be compared since they were not published. 

5.1. The comparison between Te{0~^^)opt and Te(0++)/ij 

In the presence of temperature variations across a low-density, homogeneous nebula, 
it is well known that Te(0++)opt preferentially weighs regions of high Tg, while T'e(0++)7ij 



-17- 



preferentially weighs regions of lower Tg. In such a case T'e(0++)opt > Te(0++)7ij, with 
the equahty holding in the constant-temperature case. When density fluctuations are also 
present, this basic scenario complicates, because the temperature is lower in the densest 
zones, and because coUisional de-excitations of the infrared lines may play a role. Roughly, 
wc expect higher Te{0~^~^)iR values due to the collisional suppression of the 52 /i line; as 
a result, [Te{0~^~^)opt - Te{0~^^)iFi] may become negative. Since density fluctuations are 
certainly present in the nebula, as revealed by the filaments and thin shells seen in the Ha 
images of the nebula, the value we determined for Te(0++)7R is totally plausible. 

Summing up, our Te{0'^'^)opt docs only marginally agree with the published values; 
Te{0^^)iR agrees with the most recent of the only two determinations we know of; and 
re(0++)7R compares well with the expected value of Tf,{0~^~^)opt (higher than ours by ~ 400 
- 800 K), provided there are high density zones in the nebula, which depress the fiux in 
52 II. In the following, we will list and discuss the uncertainties affecting our results. In this 
discussion we will consider several uncertainties affecting all the data, although, in the light 
of the preceding discussion and the previous comparison to published data (see Section 5), 
it seems probable that the problem lies with the adopted WFPC2 A 4363 flux, and that our 
ISO infrared data and the WFPC2 A 5007 data are reasonably accurate. 

6. Possible sources of error 

6.1. Continuum subtraction in F437N 

The long slit on which the continuum level was measured samples only one fraction of 
the whole nebula, and we cannot guarantee that such fraction is indeed representative of 
the rest. It seems not probable, however, that the continuum varies so dramatically in the 
area outside the slit, so as to alter sensibly our conclusion. A further possibility is a bias in 
the continuum level of the slit spectra. A bias in the continuum level could be caused by 
at least two different circumstances: an inaccurate correction for light dispersion inside the 
spectrograph, and scattered light from the central star within the nebula. The first effect 
could be important in view of the fact that the fiux in the continuum outweighs the fiux in 
A 4363 in the filter, so that a small bias in the continuum would be amplified in the derived 
A 4363 value. The second effect, which would mimic a more intense nebular continuum, 
would affect equally the slit spectra and the HST images, so it cannot be invoked as a cause 
for the low Te{0~^^)opt found, and will not be further discussed. 

O'Dcll & Doi (1999) stress the importance of an accurate calibration of the WFPC2 
filters in order to make a quantitative use of the HST emission-line images. One of the 
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problems they discuss related to calibration is that the filter characteristics described in 
the WFPC2 Handbook have been determined in parallel light, whereas the incident beam 
inside the WFPC2 is convergent. The consequence is that the filters as used by the WFPC2 
have shorter peak transmissions and wider FWHMs. The use of the pre-launch filters' 
characteristics provided with the HST documentation can therefore introduce a bias, and 
they recommend instead the use of calibration relations based on the instrument "as used" . 
A specific problem they address is the continuum subtraction; the strategy they suggest to 
determine the continuum level is to use an image taken in F547M, a filter that includes 
no prominent emission lines, together with an assumption on the color of the continuum. 
Unfortunately, NGC 6543 has not been observed in this filter, so this method is not applicable 
in our case. 

A different strategy to subtract the continuum from the F437N filter has been proposed 
by O'Dell, Peimbert, and Peimbert (2003) for the case of low-excitation objects. The basic 
idea is to take an image in the F469N filter: the constraint on the excitation degree of the 
object ensures that the contribution of the He II line A 4686 is negligible, and all the signal in 
the image can be interpreted in terms of continuum. This procedure allows to reduce the un- 
certainty related to the assumption on the color of the continuum, since the distance between 
the wavelengths of F437N and F469N is quite small. Unfortunately, although the He"'"+/H"'" 
ratio in NGC 6543 is very small and make this object suitable to apply the method by O'Dell, 
Peimbert, and Peimbert (2003), no F469N image of this object has been obtained (maybe it 
is precisely because the He'*"+/H+ ratio in NGC 6543 is so small, that an observation of the 
nebula with this filter was judged uninteresting by past observers). 



6.1.1. Other issues related to the HST/WFPC2 calibration 

There are further possible explanations for our anomalous Te.{0'^^)opt result related to 
the WFPC2 calibration. We will give a broad overview of them, without estimating them 
quantitatively. 

One problem is the accuracy in the determination of the contribution of the various 
lines in each filter. Due to the nature of the calibration procedure mentioned in the previous 
section, together with the natural changes suffered by any filter with time, the specifications 
reported in the WFPC2 Handbook do not necessarily give an accurate description of the 
actual filters' transmission curves. As a result, the decontamination of the F656N flux 
from the [N ll] A 6584 contribution could be slightly biased; and the decontamination of the 
F437N flux from the H7 contribution could be even more biased. Unfortunately, the relations 
provided by O'Dell & Doi (1999) to deal with this problem are not applicable, because they 
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make use of the flux in the F547M filter, which has not been used to observe NGC 6543. 

There are also other potential problems that affect only faint images with very low DNs. 
The composite F437N image used in the study has a mean DN of 2.9, corresponding on 
average to DN ~ 1.5 in each of the two exposures (the DN cannot be measured individually 
on these two images, since it is spuriously increased by cosmic rays, which are removed in the 
composite image); since these DNs arc extremely low, the effects mentioned in the following 
are possibly a large source of uncertainty in the F437N flux. On the other hand, the images 
taken with the other filters have DNs of the order of hundreds, and their calibration is much 
more accurate. 

Further relevant sources of uncertainty, which affect all CCDs at low DNs, are the read 
noise and the loss of information implied by the digitization of the signal. Another important 
source of uncertainty in our specific case is the progressive alteration in the performance of 
the WFPC2 CCDs: in orbit, the CCDs are subject to radiation damage (Biretta et al. 2001), 
which alters their long term properties leading to an increase in dark counts and a decrease 
in the charge transfer efficiency (CTE). Dark counts increase because of hot pixels, which 
are routinely identified and can, in principle, be corrected. Conversely, the decrease in CTE 
is difficult to estimate, since the existing cahbrations (Whitmore, Heyer, & Casertano 1999) 
are based on measurements of the CTE loss for a bright point source (DN ^ 60 for gain = 
7) observed against a faint background (DNjc, < 1 for gain = 7); therefore, it is not clear 
whether the correction formulae provided by Whitmore, Heyer, & Casertano (1999) can be 
equally applied to the case of a diffuse, nebular object, which can be seen as a moderately 
faint object seen against a moderately faint background (in our case, DN = DN(,g ~ 1.5 
for gain = 7). On the other hand. Figure 4.14 in Biretta et al. (2001), based on recent 
observations, includes a few more points for brighter backgrounds, in addition to those on 
which the calibration by Whitmore, Heyer, & Casertano (1999) was based. The behavior of 
the recently added points suggests that the CTE does not depend on the DN for backgrounds 
of order ^ 1.4 DN (equivalent to a DN = 2.8 at a gain of 7). Provided the CTE behavior 
is smooth with respect to time, we can roughly estimate that the CTE loss for the F437N 
filter was < 10% at the time our images were taken. 



6.2. ISO/LWS calibration 

As explained earlier, the uncertainty quoted on our adopted infrared line fluxes only 
represents the dispersion in our data sample. Additional uncertainties might come from 
the calibration itself. According to Swinyard et al. (1998), the uncertainty in the value 
of the flux for a line is given by the sum (in quadrature) of three terms: one of them 
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is the dispersion in the measured hne flux, and the other two are related to the detectors' 
responsivity cahbration. Although these terms depend on the specific observation considered, 
typical values are ~ 0.05 and 0.01. Adding them to the dispersion of our data points, we 
obtain an overall uncertainty of roughly 13% for cither detector. It is possible, however, that 
this value still underestimates the real uncertainty. Reasons for this are the following: 

1) Both the drift of the SW2 detector with time and the difference between the SWl and 
SW2 detectors are of the order of 10% of the measured flux. This uncertainty is not 
included in the previous estimation, and should be added to the overall error budget. 

2) The spectrograph has been calibrated assuming that the source is point-like and cen- 
tered on the optical axis. Since our object violates at least the first of these assump- 
tions, a correction to the flux should be done. As we stated before, such correction is 
in principle very small (of order ~ 1%). However, the beam profile of the instrument 
is not completely understood, being, e.g., much narrower than predicted by optical 
theory for reasons still unknown (Swinyard et al. 1998; Gry et al. 2001); so we cannot 
completely exclude that the extension of the source do play a role in our analysis. 

3) Most density determinations of NGC 6543 are well in excess of our value A^e — 1650 
cm-3; for example, DLW found A^e[0 ill] = 10000 cm-^ and A^e[0 ii] = 4000 cm-^ 
and Liu et al. (2001) found Ne[C\ ill] = 4800 cm-=^ and A^e[Ar iv] = 3900 cm-=^ (but, 
on the other hand, DHEW found A'e[0 ill] = 2000 cm~^). Higher density values would 
be obtained by increasing the flux of 52 fj, or decreasing the flux in 88 /i, leading in both 
cases to lower Te{0^~^)iR. 

As an illustrative example of the uncertainties affecting the ISO/LWS data, wc repeated 
the computations adopting the mean 52 /j, flux of the SWl detector instead of that of the 
SW2 detector. This different choice changes A 5007/52 by -0.05 dex and 52 /i/88 by +0.05 
dex, yielding = 8100 K and A^e = 2000 cm-^ 

6.3. Reddening correction 

As stated earlier in this paper, wc derive a C(H/9) value much lower than some of the 
values found in the literature. However, we see no way of increasing our result. Furthermore, 
even if this were the case, it would actually imply just a moderate improvement in the 
Te(0''"+)opt result. To illustrate quantitatively this point, we arbitrarily increased C(H/3) 
from 0.07 to 0.30, the largest of the values compiled by Hyung et al. (2001). Given the 
differential extinction affecting A 4363 and A 5007, this change translates into a moderate 
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decrease in the intrinsic A 5007/ A 4363 ratio, implying an increase of less than 200 K. On the 
other hand, such increase in C(H/9) would enormously affect the A 5007/52// ratio, yielding 
an increase of as much as 1600 K. 



7. Summary and conclusions 

This work was motivated by the desire to determine self-consistently the physical condi- 
tions in the bright core of NGC 6543. We did this by means of two different methods, both 
based on oxygen line ratios measured across the whole bright halo of the nebula. But, as 
the work progressed, it became evident that the archival data we used were not optimized 
for this particular task. 

The first method used is the standard nebular-to- auroral temperature diagnostic, and 
the second is the diagnostic diagram based on infrared lines developed by Dinerstein (1983) 
and DLW. The and values derived by means of the diagnostic diagram are not com- 
patible with those by DLW, but are in very good agreement to the more recent result by 
DHEW. On the other hand, the nebular-to-auroral temperature we derive is somewhat lower 
than the values published in the hterature. We discuss some possible causes for these results, 
related to the accuracy of HST, ISO, and ground-based data. 

As a sidcproduct of this work, we find that the continuum level differs from the one 
expected at Tg ~ Te(0++)opt when Lya conversion is neglected. We investigate the possibility 
that the extra continuum could arise from a high-temperature zone, possibly a shocked 
region, as suggested by the presence of X rays; however, we could not reproduce the required 
extra continuum by means of simple plasma models of the X-ray emitting bubble. An 
alternative explanation we propose is that it is generated by dust scattering of the central 
star's light. A third possibility, which is the most plausible, is that the additional continuum 
is enhanced 2-photon emission, originated by conversion of scattered Lya photons. 

In neither case is it possible to work out a theoretical prescription to compute accurately 
the continuum intensity that would eliminate the necessity of relying on spectroscopic infor- 
mation. Furthermore, either of these mechanisms would affect equally the slit spectra and 
the image, so neither would explain the marginally low temperature value that we deduce 
by means of the HST data. As a result, the most important source of uncertainty in the 
determination of the optical temperature is the continuum subtraction in the A 4363 image; 
a specific conclusion we draw is that it is not possible with these data to obtain an accurate 
bidimensional temperature map of the nebula. More generally, we review several possible 
causes of error that could be introducing a bias in our work. The calibration uncertainty of 
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WFPC2 at very low DNs introduces perhaps an additional uncertainty. The HST fluxes at 
other wavebands and the ISO data are probably more robust, but the overall uncertainty 
could still be somewhat higher than 10%. In the case of HST data, a fraction of this uncer- 
tainty might arise from the contribution of different lines to the observed fluxes; in the case 
of ISO data, the overall accuracy is probably of order ~ 13%. 

Although NGC 6543 has been previously studied by many authors, most measurements 
have been made with small apertures which sample only a small fraction of the nebula. One 
disadvantage of this approach is that it limits the possibility to compare among different 
results. On the contrary, the data wc use sample the whole region, allowing us to do a 
self-consistent analysis which eliminates the need of assumptions about the ionization struc- 
ture. Furthermore, the use of lines of the same ion allows a self-consistent determination 
of and Tg, independent of assumptions on the ionization structure. Unfortunately, the 
available archival data were not specifically designed for this analysis, and contain sources 
of uncertainties that are, in many cases, quite difficult to estimate. We think therefore that 
our results on the physical conditions in NGC 6543 are somewhat biased, and we definitely 
exclude the possibility of using these data to derive sensible temperature maps unless specif- 
ically designed data are obtained that provide a safe way of measuring the continuum and 
line contribution inside each filter. Our results are, however, largely informative for what 
concerns delicate technical issues that must be taken into account when studies of this kind 
are carried out. In particular, we draw the attention, throughout this paper, to complex 
aspects of the use of data from large archives, in the hope that the authors of future sim- 
ilar works will be motivated to provide an explicit description of the crucial assumptions 
underlying their analyses. 
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I(88/x) (10"^° erg s"^ cm"^) 

Fig. 1. — ISO 52^ and 88/^ line vs. continuum intensities for our sample, with the l-cr ellipses 
overplotted. Upper panel: dots are SWl data, open circles SW2 data. Lower panel: dots 
are SW5 data, open circles LWl data. 
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Fig. 2. — (top) Nebular spectrum of NGC 6543 along position angle 5° in the wavelength 
region around A 4363. The dashed line is the actual spectrum; the full line is the spectrum 
convolved with the F437N filter response (the total system peak response of 0.03 has been 
divided out). The inset shows a sample of the spectrum along the slit in a logarithmic scale, 
(bottom) The dotted circles (left handsidc scale) show the percent of the total intensity of 
radiation through the HST filter F437N that is contributed by the A 4363 line intensity. Also 
plotted is the variation of the A 4363 intensity along the slit (full fine, right hand side scale) . 
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Fig. 3.— Median value of [O ill] A 5007 on the HST image. 
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Fig. 4. — Hydrogen and helium continua at Tg = 8000 K. 
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QT 01 Q 

Fig. 5. — Observed H7 and nearby continuum fluxes measured along the INT slit, as a 
function of the distance from the central star. The shaded region corresponds to the region 
predicted by recombination theory. Dots and open circles represent the two directions along 
the slit, and the two solid lines the corresponding best-fit linear relations. 
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Fig. 6.— Pixel-by-pixel F(Hq;) /-F(H/3) histogram. 
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Fig. 7. — [O III] infrared-line diagnostic diagram. 
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Table 1. ISO LWS archive spectra. 



Epoch Date TDT Proposal ID 



Epoch I 



Epoch II 



05- Jun-1996 
07-Jun-1996 
14-Jun-1996 
21-Jun-1996 

27- Jun-1996 
04-Jul-1996 
ll-Jul-1996 
18-Jul-1996 
24-Jtil-1996 

28- Jul-1996 
31-Jul-1996 

02-Dec-1997 
09-Dec-1997 
16-Dec-1997 
23-Dec-1997 
30-Dec-1997 

06- Jan-1998 
13-Jan-1998 
20-Jan-1998 
27-Jan-1998 



20101905 
20301905 
21003105 
21704705 
22301705 
23001905 
23705405 
24404805 
25100805 
25500701 
25800905 

74801614 
75500604 
76201203 
76902404 
77602002 
78300202 
79000702 
79700402 
80401402 



CAL_R196 

CAL_R203 

CAL_R210 

CAL_R217 

LWS_R223 

CAL_R230 

CAL_R237 

CAL_R244 

CAL_R251 

DUST4_A 

CAL_R258 

CAL_R748 
CAL_R755 
CAL_R762 
CAL_R769 
CAL_R776 
CAL_R783 
CAL_R790 
CAL_R797 
CAL_R804 
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Table 2. ISO/LWS fluxes. 



Line intensity^ Continuum fiux^ 



Line 


Detector 


Epoch^ 


Median 


a 


Median 


cr 




SWl 


I 


5.58 


0.50 


1.66 


0.29 






II 


6.10 


0.43 


1.91 


0.68 






I+II 


5.74 


0.53 


1.70 


0.50 




SW2 


I 


5.06 


0.66 


1.54 


0.11 






II 


5.15 


0.52 


1.49 


0.12 






I+II 


5.10 


0.61 


1.53 


0.12 


88/i 


SW5 


I 


1.40 


0.18 


2.69 


0.47 






II 


1.42 


0.13 


2.89 


0.29 






I+II 


1.41 


0.16 


2.79 


0.41 




LWl 


I 


1.34 


0.14 


2.92 


0.19 






II 


1.44 


0.20 


3.09 


0.17 






I+II 


1.39 


0.17 


3.02 


0.20 



^In units 10 ^° erg sec ^ cm ^. 

''In units 10"^° erg sec"^ cm-^ ^-i for SWl and SW2, and 
10-^1 erg sec-^ cm'^ //-i for SW5 and LWl. 

^See text and Table 1 
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Table 3. HST archive images. 



Line^ 


Filter 


Dataset 


Exposure 


PHOTFLAM 


RECTW 








(s) 


(erg cm~^ A~^) 


(A) 


[0 III] A 4363 


F437N 


U27Q0103T 


1200 


7.400 X 10-16 


31.829 






U27Q0104T 


1200 






H/3 


F487N 


U27Q0105T 


500 


3.945 X 10-16 


33.921 






U27Q0106T 


900 






[0 III] A 5007 


F502N 


U27Q0107T 


200 


3.005 X 10-16 


35.781 






U27Q0108T 


600 










U27Q0109T 


200 










U27Q010AT 


600 






Ha 


F656N 


U27Q010FT 


200 


1.461 X 10-16 


28.336 






U27Q010GT 


600 






[N II] A 6584 


F658N 


U27Q010HT 


500 


1.060 X 10-16 


39.232 






U27Q010IT 


500 







^Most important emission line included in the passband. 
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Table 4. Frequency dependence of some continuum-emission coefficients (in 10 erg 

sec"^ cm^ Hz-i) at = 8000 K. 





A 4363 


A 5007 


A 6563 


H P 


2.196 


3.687 


8.431 


He P 


2.544 


4.502 


10.12 


2-pfioton continuum, minimum'' 


3.827 


3.193 


2.175 


2-pfioton continuum, maximum^ 


13.429 


11.204 


7.632 



^Including recombination and bremsstrahlung. 
^Neglecting LyQ;-to-2-photon conversion. 
^Including maximal LyQ;-to-2-photon conversion. 
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Table 5. Summary of fluxes and diagnostics. 



DLW This work 





F(A) /(A) 


F(A) /(A) 


A 5007^ 


1.65 ±0.25^ 2.70 ±0.40 


1.65^ 1.93 ±0.29 


A 4363^ 




3.28 ±0.40 3.95 ±0.84 




8.79 ± 0.97 


5.10 ±0.61 




1.10 ±0.20 


1.39 ±0.17 


log A 5007/A 4363 


2.56 


2.69 ±0.10 


log A 5007/52/x 


49+0-07 




log 52/i/88/x 


90+0 09 


S6+0-07 




7900 


7500 ± 450 


Te(0++),« 


5800 ± 300 


8600 ± 500 


log N,{0++),R 


4 0+°° 


3.22 ±0.12 



^In units of 10 ^ erg s ^ cm ^. 

^Computed by us with the C(H/3) quoted by DLW. 

■^No uncertainty have been quoted on F(A5007), since our estimated 
uncertainty on F(A5007) includes only the uncertainty on the reddening 
correction (Section 5). 

'^In units of 10~^^ erg s~^ cm~^. 



